r 



I 

i t 

i 

i 

1 


\ 


^Ml- ANNUAL REPORT 
Grant No. NAG-849-1 


p </? 


PREDICTION OF FORCES AND MOMENTS 
FOR HYPERSONIC FLIGHT VEHICLE 
CONTROL EFFECTORS 



-i 


Mark D. Maughmer 

Associate Professor of Aerospace Engineering 


Lyle N. Long 

Assistant Professor of Aerospace Engineering 


Peter J. Pagano 
Graduate Research Assistant 


Department of Aerospace Engineering 
233 Hammond Building 
The Pennsylvania State University 
University Park, PA 16802 
(814) 865-2569 


submitted to 


National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 
Attn. David Raney 


(f.ASA-C K - 1 ^ fl 9 *> A ) PREDICTION Of FH^Cf.S AND N92-11070 

MOMENTS F^P hYPt RSClNI C FLIGHT VTliICL' 

COMROL EFFTCTORS Semi annual Report 

(P^nnsvl vani a State Univ.) ab p CSCL 223 Unci as 

G3/1H 00A6A6A 


I 

I 

t 



TABLE OF CONTENTS 




SUMMARY 

1. INTRODUCTION 

Purpose of this ./ork 
Methodology 

2. THEORETICAL ASPECTS 

TEAM Code 

Surface Inclination Methods 

Modified Newtonian Theory 
Tangent Wedge Method 
Hypersonic Shielding 

3. TEAM CODE VALIDATION 

Supersonic/Hypersonic Inviscid Cases 
Viscous Cases 

4. PROGRESS 

Phase I 
Phase II 

5. CURRENT AND FUTURE WORK 

6. REFERENCES 

APPENDIX. TEAM Code Description 
FIGURES 


PAGE 

2 

2 

2 

3 

4 

4 

5 

5 

6 
7 

S 

S 

9 

10 
10 
10 

15 

16 
20 
25 




SUMMARY 

The first phase of this research effort, developing methods of predicting flight control forces and 
moments for hypersonic vehicles, included a preliminary assessment of subsonic/supersonic panel 
methods and hypersonic local flow inclination methods for such predictions. While these findings 
clearly indicate the usefulness of such methods for conceptual design activities, deficiencies exist 
in some areas. Thus, a second phase of research was proposed in which a better understanding is 
sought for the reasons of the successes and failures of the methods considered, particularly for the 
cases at hypersonic Mach numbers. This second phase was begun in March, 1990. 

To obtain this additional understanding, a more careful study of the results obtained relative to the 
methods utilized was undertaken. In addition, where appropriate and necessary, a more complete 
modeling of the flow has been performed using well-proven methods of computational fluid 
dynamics. As a result of this more complete understanding, assessments will be made which are 
more quantitative than those of Phase 1 regarding the uncertainty involved in the prediction of the 
aerodynamic derivatives. In addition, with improved understanding, it is anticipated that 
improvements resulting in better accuracy will be made to the simple force and moment prediction 
methods considered in this study. 

1. INTRODUCTION 

Purpose of This Work 

Included in the executive summary of the AGARD Symposium on the Aerodynamic Characteristics 
of Controls in 1979 [1] was the need for a more extensive and modern data base. Furthermore, it 
was suggested that additional research be conducted to fill gaps in the data base. It was also 
pointed out that theoretical methods were inadequate in accounting for viscous effects and flow 
separation. More than ten years later, these comments still apply. 

In fact, for hypersonic flight vehicles the situation is actually worse. The data base for hypersonic 
flight control information is extremely limited. Some available wind-tunnel data is of questionable 
validity and flight-test results are scarce. Furthermore, in addition to the need to account for 
viscous effects and flow separation, theoretical prediction methods at hypersonic Mach numbers 
must also contend with problems involving thin shock layers and real gas effects. Some existing 
computational fluid dynamic (CFD) methods have the ability to handle such problems, but require 
too much computer and engineering time to be used routinely for conceptual design studies. 
Consequently, a real need exists for computationally efficient methods of predicting flight control 
forces and moments for hypersonic vehicles which nevertheless provide reasonable results. 

The recent push toward the development of hypersonic flight vehicles has highlighted the need for 
rapid aerodynamic prediction methods (2 - 4], The complexity of such vehicles demands the 
integration of all technological disciplines from the conceptual design stage. Hence, it is of great 
advantage to be able to analyze many conceptual design proposals and discard those which are not 
promising in a timely manner. Many methods exist which are capable of performing this analysis, 
however, at the conceptual design stage monetary and/or time restrictions may preclude their use. 
CFD techniques are best suited for preliminary or detailed design analysis due to the great length 
of time required for solution of flows over complex geometries. In addition, the expense and 
limitations of hypersonic test facilities may relegate their effective use to the testing of final design 
configurations. 
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Early integration of control systems into the design process is of paramount importance to the 
success of a hypersonic vehicle design. At hypersonic Mach numbers, a vehicle traveling ‘trough 
the upper atmosphere will experience dissociation of constituent gases in air. The Space Shuttle 
Orbiter is a case in point. Upon re-entry, STS-1 required a body flap deflection twice that of the 
predicted value to trim out the longitudinal moment [5 - 10). While there is some disagreement 
over the cause of this problem, most believe it to be due to either real gas effects [11, 12] or low 
Reynolds number effects [13]. This question is being examined more closely. Other complicating 
features of hypersonic flows are thick boundary layers, entropy layers, thin shock layers, and 
boundary layer/shock layer interaction, all of which effect the control aerodynamics. Hypersonic 
vehicles also experience very large center-of-pressure movements as they traverse the flight 
envelope from low to high speed. Likewise, the design engineer must consider the changes in nap 
effectiveness due to the flap being embedded in the boundary layer. It is clear from entropy layer 
studies [14] that a sharp nosed cone produces a greater pressure recovery on a deflected flap surface 
than a blunt one. Thus it is concluded that flap effectiveness is decreased by increasing nose 
bluntness. It is interesting to note that widely used aerodynamic prediction techniques, such as the 
Aerodynamic Preliminary Analysis System (APAS), other Gentry codes and the Hypersonic 
Arbitrary Body Program (HABP) make no attempt to model the flow field ahead of a control 
surface. Flap effectiveness is also decreased as the flap deflection angle is increased. This is 
caused by the boundary layer separating when the flow is deflected far enough relative to the 
vehicle body. This separation can create a secondary shock system and/or transition the boundary 
layer, both of which decrease flap effectiveness. It has also been pointed out that there is a need 
for a substantial data base [14], which must be shared among the various disciplines involved in 
hypersonic research. 

The purpose of this work is to: 1) to develop and/or improve pie hypersonic aerodynamic 

methods such as those used in APAS/HABP to better predict viscous and dissociating flow field 
effects on control surfaces. The initial phase [15] consisted of comparing these simple prediction 
techniques (data generated using APAS) to experimental data to determine where the greatest need 
for improvement lies. In this second phase we have been trying »o determine which flow 
phenomena (ie. shockwaves, viscosity, chemistry) have the greatest effect on prediction quality. 
2) Add to the database of hypersonic research by examining a relatively simple geometrical shape 
(the X-15 airfoil) for a range of Mach number, angle of attack, flap deflection and flow field 
conditions. And 3) attempt to clarify the "hypersonic anomaly" experienced by STS-1. 

Methodology 

Portions of the three goals of this research project, as outlined in the preceding section, are being 
attained through the use of an advanced CFD code (the TEAM code), see Chapter 2. This code 
is essentially being used as a hypersonic wind tunnel. CFD is playing a very important role in the 
advancement of hypersonic research. Flight testing is a valuable means of collecting data, but it 
is difficult to accomplish and is often performed post-development Experimental ground facilities 
are simply too limited to cover the range of parameters and flight conditions [16]. CFD is not 
currently the complete solution, there are still many problems to be overcome. It has been 
suggested, however, that it may provide results as meaningful as those obtained from experimental 
ground facilities. Such facilities are plagued by the need to extrapolate data to flight conditions, 
contaminated flow [5], and tunnel peculiar effects on produced data [14). A major contributor to 
the tunnel peculiar effects is that of the acoustic environment. The active turbulent boundary layer 
on the wall of a hypersonic tunnel, as well as any other acoustic disturbances of sufficient strength 
introduced into the flow-field, will cause transition to occur on the model earlier (at a lower unit 
Reynolds number) than would be the case for a free-flight experiment [12 - 21] 
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The lack of experimental hypersonic facilities is yet another impetus for the development and use 
of CFD codes A 1968 report to the NASA Subcommittee on Fluid Mechanics of the Committee 
on Basic Research [22] mentions the need for "wind-tunnel facilities with higher Reynolds number 
capabilities than are currently available". All the more poignant in the 1990 s after decades of 
inactivity in hypersonic research have depleted the number of operational test facilities The high 
cost of hypersonic test facilities has often been their demise. Reference [23] quotes the costofa 
Re (based on test section area) * 10 million continuous flow tunnel as over 100 million 1975 dollars. 

Thus with an understanding of the shortcomings associated with CFD, and an appreciation of its 
advantages, a comparison is presented among results from classical Newtonian and tangent wedge 
theories (both integral components of the APAS/HABP type codes) and results obtained from the 
TEAM code. Surface pressure plots are compared to illuminate anomalies, and flow-field contour 
plots from the TEAM code will be shown to explain these differences. Then recommendations will 
be made for the improvement of the APAS/HABP type analysis methods. 

2. THEORETICAL ASPECTS 

As discussed in Chapter 1, a comparison is presented between simple hypersonic methods for 
predicting surface pressure and an advanced CFD technique. The theory pertaining to these various 
methods is presented here. 


The Three-Dimensional Euler/Navier-Stokes Aerodynamic Method (TEAM) [24] was developed by 
the Lockheed Aeronautical Systems Company, Burbank, California for the Aeromechanics Division 
of the Flight Dynamics Laboratory, Wright Research & Development Center under contract of the 
United States Air Force. The four year (July 1984 - October 1988) effort is the result of a desire 
to develop a computationally efficient code which could solve both viscous and »nviscid How fields 
with real-gas effects. Grid system independence was another driving factor in TEAM code 
development. Grids may be generated using any external program available to the user, only the 
cartesian coordinates of the grid nodal points are required. 

The TEAM code will be briefly describe here, a more detailed description is included in the 
Appendix of this report. TEAM uses a finite-volume spatial-discretization algorithm coupled to 
a Runge-Kutta time-marching scheme to solve the Navier-Stokes equations. It can use zonal, 
patched grids and is therefore quite flexible. It also includes implicit residual smoothing, enthalpy 
damping, and local time-stepping for efficient convergence to a steady state. The code can 
simulate turbulent, laminar, or in viscid flow of perfect or real gases. 

The TEAM code is normally run on a large supercomputer such as a Cray, Convex, or IBM-3090. 
As a rough rule of thumb, the code requires 25 microseconds/grid point/iteration. Just to illustrate 
the CPU requirements of this code, some representative cases are illustrated in the following table. 
This is just an illustration and individual times for particular runs can vary dramatically depending 
on the flow field and the grid. 
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Dimensions 


No. of Cells 


No. Time Steps 


CPU Time 
(Cray XMP) 


[ Inviscid (Euler) 

2-D 

3-D 

10,000 

400,000 

1,000 

2.000 

4 Min. 

5 Hours 


Viscous (Navier-Stokes) 


3-D 


20.000 


10.000 


1.4 Hours 


800,000 


20.000 


111 Hours 


This table shows that 3-D, viscous computations are extremely time consuming. It also shows, 
however, that even a 2-D viscous computation can require hours on a supercomputer. This is one 
reason for performing the work described herein. In a design environment, one must have tools 
which are computationally efficient in order to facilitate the iterative nature of design. 


The nonlinear nature of hypersonic flow manifests itself in such phenomena as high-temperature 
chemically reacting flow-fields, thin shock layers, entropy layers (vorticity interactions), 
interactions between the viscous boundary-layer and the shock wave and low-density effects at high 
altitudes. Considering these phenomena, to obtain a complete picture of the flowfield, one cannot 
hope to use a simple analytic method. However, there do exist a number of analytic methods, 
which under certain circumstances provide a good first approximation to the coefficient of pressure 
(and hence the aerodynamic forces and moments) on a body in a hypersonic inviscid flow-field. 

Modified Newtonian Theory 

Sir Isaac Newton developed his famous Newtonian flow model more than three centuries ago, it was 
first published in Propositions 34 and 35 of Principia in 1687 [32] Although developed to explain 
subsonic flow, this method has seen renewed interest in the latter half of this century as a means 
of predicting the aerodynamic forces on hypersonic vehicles in the design proposal stage. The 
equation for this model may be obtained by assuming that a flow impacting on a surface would 
loose all momentum normal to the surface, and the flow particles would then move tangentially 
along the surface. Then for a surface inclined at an angle 6 to the free-stream. 

Change in normal velocity = V m SIN6 
Mass flux incident on a surface A = p.V.A SIN6 

Time rate of change of momentum of the mass flux - p m K,ASIN 0 

Newton’s Second Law states that the time rate of change of momentum is equal to the force exerted 
on a surface, denoting this force by F, 

F = p.vf/4sin J 6 


Newton assumed the flow of particles to be rectilinear, ie. no random interaction of fluid particles, 
hence F is associated only with the linear motion of the particles. Static pressure of a gas is due 
to the purely random motion of its particles, not accounted for in the Newtonian model. Thus F/A, 
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free-stream static pressure. 
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This result can also be derived from .he gas dynamic equations governing oblique shock waves. 

normal shock. 

. In (2.15) 
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Newtonian in the prediction of pressures over blunt bodies. 


Tangent Wedge Method 

The .angem wedge method was developed ,o predict surface The 

dimensional hypersonic shapes. Referring to ig • * surface of the bodv at point 
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to the free-stream be greater than the maximum 
The tangent wedge method has been shown to 
leading edge shocks. 


turning angle for the free-stream mach number, 
work best for sharp nosed bodies with attached 
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Hypersonic Shielding 

Hypersonic shielding is a method used to treat the leeward side of bodies in 

These wake regions are "shielded" from the oncoming free-stream, as a result the surface 

mav be set to zero static pressure. Hypersonic shielding is often used as it is in this research 

effort, in conjunction with Newtonian theory, or the tangent wedge method for the impact side of 

the body. 


The definition of pressure coefficient is 


so when p - 0 (vacuum) : 
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and since * RT » — , the above becomes : 
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which is the pressure coefficient in a vacuum. 



3. TEAM CODE VALIDATION 

Experimental verification of the numerical results obtained in this research has not been possible 
due to an absence of access to a hypersonic wind tunnel facility. Validation of the TEAM code 
for similar (ie. hypersonic in viscid and viscous) experiments will therefore be cited as verification 
of the TEAM code. 


Sypersonlc/Hvpersonic Inviscid Cases 

The first test case to be discussed involves evaluation of TEAM’S ability to model attached and 
detached shocks in supersonic flows. This was accomplished by running TEAM for both a sharp- 
nose cone cylinder and a blunt-nose cone cylinder at a Mach number of 2.96 and angles of attack 
of zero and 16 degrees using the standard adaptive dissipation scheme. Additionally, the sharp- 
nose cone cylinder was modeled at Mach number 4.63 for angles of attack of 4 and 24 degrees. 
Surface pressure data obtained from TEAM is correlated with experimental data from reference 
[35]. 

Good agreement between the experimental and computational results is apparent. Computational 
and experimental data are again in good agreement except at 6 = 45°, here computed surface 
pressures are below the measured results. The authors of reference [24] cite the absence of viscosity 
as the most likely cause for the variance. 


Evaluation of TEAM’S ability to predict hypersonic flow-fields is accomplished by trying to 
duplicate the experimental results of Shindel [36]. A cone-derived hypersonic waverider was tested 
at Mach = 6 for angles of attack of -4°, the design angle of attack of 0°, and +4°. Computed lift 
and drag coefficients were in very good agreement with experimental and theoretical values [36]. 
The data on the upper surface and lower surface out to about 60% correlate well. The disagreement 
beyond this point on the lower surface is primarily due to the absence of a shock/boundary layer 
interaction in this region. 


A further investigation of hypersonic waverider configurations at off design conditions using the 
TEAM code has been performed by Long [37], Inviscid perfect and real gas computations were 
performed on Rasmussen’s elliptical-cone waverider [38]. Computed values of CL, CD and L/D 
at angle of attack correlated very well with experimental data. Differences between experimental 
and computational results were attributed to the lack of viscous effects in the computational effort. 
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Viscous Cases 

The first viscous case to be discussed is that of the Lockheed-AFOSR Wing C. This test was 
conducted at Mach = .85, a * 5° and R e * 10 million based on a MAC of .76 meters. The thin- 
layer (TLNS) approximation to the RANS was solved in the six zones surrounding the wing, and 
the Euler equations in the seventh outer zone. Computed surface pressure data are compared with 
experimental results [39] and an inviscid solution at four stations along the wing. It is evident that 
the viscous solution better approximates the experimental data. The viscous solution shows a 
reduction of aft loading and a forward movement of shocks on the wing. It is stated in reference 
[24] that further research is required to explain the differences between the viscous solution and 
the experimental data, especially at the outboard stations. 


A viscous solution for a double-delta wing- body configuration is another case that was analyzed 
for the TEAM code validation. This case was run at Mach = .3, a = 20° and R e = 1 million/ft. 
The TLNS equations were solved in 12 zones next to the body, the Euler equations in the remaining 
four. Turbulent flow was assumed over the entire surface. Surface pressure data at three cross 
plane stations for TEAM viscous, TEAM inviscid and experimental are presented in Reference 38. 
Better over-all prediction at all stations is evident for the viscous solution. Lift and drag 
coefficients for the viscous case also better approximate the experimental data [24], 

Reference [40] is a viscous investigation of an axisymmetric indented nose cone at Mach = 9.89, 
intended to compare the results from continuum (Navier-Stokes, Euler) methods with those of a 
kinetic theory approach (Boltzmann equation). The Navier-Stokes method used is TEAM. The 
kinetic theory approach, as mentioned above, is the DSMC (Direct Simulation Monte Carlo) method 
developed by Bird [41], Heat transfer predictions for each method are shown to compare very well 
to experimental data from reference [42], Other comparisons of surface pressure coefficient, skin 
friction coefficient, flow-field density and temperature are also made between the two methods 
without comparison to experimental data. These flow-field correlations are quite good, and show 
at least a good agreement between the methods used in TEAM and a solution of the Boltzmann 
equation. 
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4. PROGRESS 


f h nt I 

To help address the need for flight control prediction tools, a research program has been underway 
specifically to provide methods suitable for conceptual design activities involving aerodynamic 
flight controls. The initial phase of this research included cataloging existing data for hypersonic 
vehicles and comparing these data with computationally efficient prediction methods. In particular, 
a preliminary assessment of the subsonic/supersonic panel methods and the hypersonic 
Newtonian-flow based methods incorporated in the APAS/HABP code [43, 44] has been made. 
This assessment [45] included a comparison of theoretical predictions with results obtained 
experimentally for the flight vehicles: the North American X-15, the Hypersonic Research 
Airplane, and the Space Shuttle. While the experimental data used was taken primarily from 
wind-tunnel measurements, a few flight-test results for the Shuttle were also included. 
Comparisons were made from Mach numbers of near zero to twenty. 

It was shown that the flow inclination methods do a good job of predicting lift at hypersonic 
speeds. Most important for flight controls work, the change in lift coefficient due to an elevon 
deflection is predicted very well. However, the pitching moment versus angle of attack and control 
deflection angle for the Shuttle at a Mach number of 5.0 is not predicted nearly as well as the lift 
curve predictions. The change in pitching moment coefficient with control deflection is reasonable, 
especially at the higher angles of attack. Since separation is not modelled, the results for large 
control deflections are not good at all. Predicted lateral/directional results (at a Mach number of 
5.0) not only agree well with wind-tunnel results, but also agree reasonably well with flight-test 
data. 

This Phase I work has been widely reported. This documentation includes two Master’s theses [45, 
46], two NASP Contractor Reports [47, 48], and a conference paper [49]. 


Phase II 

During Phase II, the reasons why the impact methods do not accurately predict some aerodynamic 
coefficients are being explored. The Three-Dimensional Euler/Navier-Stokes Method (TEAM) 
developed by Lockheed for the US Air Force is being used to accomplish this. This code has been 
thoroughly tested and can be used as a numerical wind tunnel. It uses a finite-volume, Runge- 
Kutta algorithm to solve either the Euler equations or the Navier-Stokes equations. It can also 
model real gas effects and turbulent flow. 

To do this, experimental forces and moment data is not sufficient to explain the differences 
between the impact methods, rather, surface pressure data is required. As experimental data such 
as this is essentially non-existent for hypersonic Mach numbers, advanced CFD techniques are being 
used to generate it. 

Initially, the X-15 airfoil is being used as a test case. This is a modified NACA 66-005 airfoil. 
For the purposes of this study, the NaCa 66-006 airfoil is scaled down to 5% thickness, since the 
NACA 66-005 coordinates are not readily available. Unlike the actual X-15 airfoil, the leading 
edge radius was not scaled and the trailing edge was not blunted Thus far. at a Mach number of 
6.93, angles of attack of 0, 10, and 20 degrees have been considered, with flap deflections of -10, 
0, and 10 degrees. At M=23, alpha « 0 and 30 degrees with no flap deflections have been 
considered. 
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Currently, results for this airfoil have been obtained from both the above described CFD code (in 
inviscid mode) and from impact theory. Since this is a 2-D problem, instead of using HABP, a 
simple program developed by Dr. Long is being used. This program, herein called HyperAero, 
simply uses the modified Newtonian method, the tangent wedge method, and a vacuum condition 
for surfaces facing away from the free-stream In the future. APAS/HABP will be used. 
However, using very simple methods initially allows the determination of what physics must be 
modelled better in order to more accurately predict the forces and moments. The HyperAero 
program simply uses the Modified Newtonian method or the Tangent Wedge method for exposed 
panels and vacuum for hidden panels. 

The grid used in the TEAM code is shown in Figure 1 and had 6,144 cells. Several grid sensitivity 
studies were conducted to determine how fine to make the grid. The grid used is extremely fine, 
but should results in accuracies within a few percent. 

W'hile the above methods are the simplest approximations possible for hypersonic aerodynamics, 
their limitations are easily quantified. Thus, a great deal can be learned by comparing modern CFD 
methods to the above methods. In addition, because HyperAero is a Fortran program that is only 
about ! 50 lines long, modifications and numerical experiments can be performed very easily. Once 
it is understood how to obtain accurate predictions for control surface deflections, these can be 
incorporated into APAS/HABP. 

Surface pressure, grey-shaded flow field images, and aerodynamic coefficients are presented below. 
Figures 2a - 12a show surface pressure predictions (from HyperAero and TEAM) for the Mach = 
6.83 inviscid cases. Figures 2b - 12b show the TEAM code flow fields shaded according to 
pressure for the Mach = 6.83 inviscid cases. The shading actually corresponds to ln(p), in order 
to show the gradients better. Since HyperAero predicts only surface quantities (not the whole flow 
field), the TEAM and HyperAero flow fields cannot be compared directly. 

Figures 2a, 3a, and 4a show the surface pressure predictions from TEAM, tangent wedge, and 
Newtonian flow for Alpha = 0 Qegrees and de!ta=-)0, 0, and +10 degrees, respectively, where delta 
is the flap deflection angle. The Euler code solution is quite smooth, but the HyperAero solution 
shows a rapid change in Cp at the mid-chord. This is due to the sudden change from the Tangent 
Wedge method to the vacuum method. In the Tangent Wedge method C - 0 when 0 = 0; 

however, as soon as 0sO the method gives the vacuum value for Cp. It is interesting to note that 
the Euler code produces a smooth result that is approximated by the HyperAero code. It should 
also be noted that the Euler code predictions are very close to vacuum conditions on the upper 
surface near the trailing edge, so HyperAero and TEAM agree quite well there. Also, the 
HyperAero method agrees well with TEAM on the upper surface of the flap since the flow is quite 
well approximated by a vacuum in that region. The Hyperaero predictions on the lower surface 
of the flap do not agree well with TEAM, however. The tangent wedge method is quite different 
than TEAM. This severely effects the moment prediction. The tangent wedge method, however, 
agrees quite well with TEAM near the leading edge. 

The impact methods are known to predict low angles of attack poorly. In fact, Hankey (501 claims 
that more refined methods must be used for angles of attack below 10 degrees. He also says that 
"only gliders having L/D’s greater than 4 will fly at angles of attack less than 10 degrees." So at 
low angles of attack, one can expect problems in using impact methods. 




Figure 5a shows the results of the three prediction methods (TEAM, tangent wedge, and 
Newtonian) for alpha=10 and delta=0 degrees In this case, the tangent wedge method agrees better 
with TEAM than the Newtonian method, which is what one would expect. Even though the 
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tangent wedge method predicts the surface pressure very well over most of the airfoil, however, f 

differences are observed between the TEAM code and tangent wedge at the leading and trailing t 

edges. This immediately indicates a possible error in the moment prediction. The leading edge ! 

discrepancy is due to the transition from tangent wedge (or Newtonian) to the vacuum method. ' 

The flow transitions from compression to expansion very near the leading edge, but the jump from < 

Newtonian to vacuum is not as noticeable at alpha* 10 degrees as it was at alpha-0 degrees. While 
the TEAM code gradually expands the flow around the corner, the impact methods transition quite 
abruptly. 

In figure 6a the results for alpha-10 degrees and delta-10 degrees are shown. While tangent wedge J 

clearly agrees better with TEAM over most of the airfoil, the Newtonian method once again agrees 

better over the flap surface. In figure 7a, the results for alpha-10 degrees and delta--I0 degrees j 

are shown. In this case, the flap is aligned with the flow direction, but the airfoil is at an angle j 

of attack. Since the impact methods do not account for any upstream influences, the flap does not 

"know" the airfoil is at alpha-10 degrees. Consequently, the tangent wedge and Newtonian methods 

both predict the same Cp’s on the flap that they predict for alpha-0 degrees and delta-0 degrees. 

The TEAM code flow field, however, has already been turned by the airfoil and must re-expand 
(lower surface) or re-compress (upper surface) the flow. This discrepancy in the impact methods 
may cause a significant error in the moment and, possibly, even the lift. -j 

Figure 8a shows the results for alpha-20 degrees and delta-0 degrees. At these higher angles of 

attack, the impact methods become quite effective and the tangent wedge method clearly agrees 

better with TEAM. However, you can still see the discrepancies at the leading edge due to the <J 

rapid transition from tangent wedge to vacuum method. There are also inconsistencies at the 

trailing edge. 
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Figure 9a and 10a show the surface pressure predictions for alpha-20 degrees and delta-10 degrees 1 

and delta-- 10 degrees, respectively. For delta- 10, neither the Newtonian or tangent wedge method i 

agree with TEAM on the flap. This case results in very strong shockwaves, one starting at the , 

leading edge and the other starting at the flap leading edge. These are clearly visible in the flow ! 

field shown in Figure 9b. It will be important to obtain viscous results for this case, since 
significant shock/boundary layer interactions probably exist. The Euler results cannot be assumed 
correct for this case. 

For the negative flap deflection (Figure 10a), the tangent wedge method does a fairly good job ’ 

predicting the TEAM code results. 

Figures 1 1 and 12 show predictions for Mach = 23. At alpha = 0 degrees, the methods all agree 
fairly well except at the leading edge. At alpha-30 degrees, the shockwave is extremely strong and 
both Newtonian and tangent wedge deviate from the TEAM code results. It will be important to 
evaluate the real gas effects at this Mach number. 


In all the flow field images, one can clearly see the shock wave starting at the leading edge. For 
a point .source travelling at Mach « 6.83, one would expect a Mach cone angle of 

P - — * — — * 0.145 radians 
M 6.83 

P * 8.4 degrees 

This is roughly the shock angle near the nose shown in Figure 2b (alpha « 0 degrees). If the airfoil 
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31 


under consideration were a flat plate, the relation between the shock angle, beta, and the angle of 
attack, alpha would be governed by 


p ? +1 * !n!2 l * 1 

« 1 T ^ 4 H* a 1 


So for alpha - 10 degrees, Mach - 6.83, and gamma - 1.4, one would get: 

P « 16.3 degrees 

So the ancle between the airfoil and the shock would be roughly 6.3 degrees. This is roughly what 
fs observed in Figure 5b. The Newtonian flow method assumes the shock wave lies directly 
body surface, which is only true for infinite Mach number and gamma* . 

The surface pressure plots and flow field images presented in this section provide an enormous 
ImounToTdau From^ the above d, session one can see how readily the source of .he errors .n the 
pnpaci' methods can be determined. During Phase 1 of .his gran, numerous problems we e 
uncovered in .he impact methods, but force and moment data alone is not enough mformatto 
to propose modifications to the impact methods 01 their usage. 

The force and moment predictions for the preced ng cases are summary i" 
is often difficult digesting large tables of numerical data, ,t is included I f« 

data illustrate how well the force and moment data can agree, even though the surface pressure ma> 
b^poo^Tpredicted. They also show how, in other cases, the different prediction methods can 

disagree by a large margin. 

In conclusion one must be careful in interpreting the results presented here. So far, only inviscid 
fi e Euler equation) TEAM code results have been presented. One cannot assume that these are 
100% correct In some of the cases presented, viscous effects may be quite significant. These will 
“SiS.. progress report and are being run a, this time It w,l "*££££ 
compare all the methods. While it is quite legitimate to compare inviscidmethodstotag 
wedge Newtonian, and vacuum techniques, since they can all be derived from "nviscid gas 
dynamics- methods such as modified Newtonian have some empiricism built into them. Thus, while 
/will be'strictlv fortuitous when methods based solely on inviscid techniques agree with viscous 
results; the empirically based methods may be expected to predict some viscous behavior since they 
have been "tuned" to yield correct results. 
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Table 1. Forces and Moment Predictions 


Newtonian 
Tang. Wdg. 
TEAM 
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5. CURRENT AND FUTURE WORK 

Thus far, APAS/HABP, TEAM, and HyperAero for hypersonic flows have been used in this study. 
The early part of Phase II consisted of setting up the programs and gaining experience in their use. 
Now that the codes are all in place, with color graphics interfaces, it is expected that the progress 
will be significant and steady. As described above, investigations using the Euler code results for 
the study of forces and moments due to control deflections are underway. In the near future, 
Navier-Stokes (i.e. viscous) and real gas effects will be explored. Unlike experiments or flight 
tests, CFD allows viscous and real gas effects to be included one at a time or simultaneously. 

Also, much of the corroborative development of experimental, analytical, and CFD that took place 
in the other speed regimes never took place for hypersonic flows because of the hiatus of work in 
hypersonic flow. There are few detailed comparisons between the simple hypersonic impact 
methods and CFD -- especially for control deflections. 

The space shuttle experienced a severe pitching moment that was not predicted before the first 
flight. This could have caused a fatal accident if large body flap deflections had not been possible. 
This problem will be investigated using TEAM, APAS/HABP, and HyperAero to determine why 
the moments were so poorly predicted. This will also demonstrate how to model the critical physics 
in APAS/HABP in order to simulate such flows. 
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APPENDIX: TEAM Code Description 


The Three-Dimensional Euier/Navier-Stokes Aerodynamic Method (TEAM) [24] was developed by 
the Lockheed Aeronautical Systems Company, Burbank, California for the Aeromechanics Division 
of the Flight Dynamics Laboratory, Wright Research & Development Center under contract of the 
United States Air Force. The four year (July 1984 - October 1988) effort is the result of a desire 
to develop a computationally efficient code which could solve both viscous and inviscid flow fields 
with real-gas effects. Grid system independence was another driving factor in TEAM code 
development. Grids may be generated using any external program available to the user, only the 
cartesian coordinates of the grid nodal points are required. Furthermore, a grid may be subdivided 
into multiple zones, each zone having its own topology, as well each zone may be speci le or 
solution by a different method, ie. "zone 1" may be solved using the Reynolds averaged^ Nav ' e '* i 
Stokes equations, and “zone 2" by using the Euler equations, etc. Zones are specified and "patched 
together with a boundary condition data file which is read by TEAM at execution. 


The Reynolds averaged Navier-Stakes equations (RANS) are the widely used mathematical models 
for the flow of a turbulent gas in thermodynamic equilibrium. These are the equations for 
conservation of mass, linear momentum and energy which have been time averaged. In integra 
form. 


IftfdQ + (f‘ dA - ff’dA. 

K 1 a 
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where w is the vector of non-dimensionalized dependant variables 
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F* and F v are the convective and 


viscous flux vectors given by, 

l *“■ \ 


F‘ = 


PV„ ♦ P n, 

pWu. 
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Here p is the mass density, u> are the three cartesian velocity components, E is the total energy, n 
denotes a unit normal vector to the surface, q*, is the non-dimensional free-stream speed and Re^ 
is the free-stream Reynolds number based on a characteristic length. A subscript n means the dot 
product with the vector n has been taken. Standard summation notation is employed with the 
subscripts i (or j or k) - 1, 2, and 3, which correspond to the cartesian coordinates X, V, and Z, 

respectively. //»£-- is the total enthalpy, and P is the static pressure. For Newtonian fluids, 

the viscous stress tensor, r^, and heat flux, qj. are related to the mean flow quantities by. 
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where T is the static temperature, Pr is the free-stream molecular Prandtl number, y is the ratio 
of specific heats, /i e and k t are the effective dynamic viscosity and thermal conductivity 
respectively, each is the sum of z molecular and a turbulent part. The effective secondary 
viscosity, A e , is defined by Stokes’ hypothesis to be -2/i e /3. Stokes’ hypothesis is most accurate for 
a monatomic gas, and is a reasonable approximation for incompressible air; however, it’s use in 
compressible applications is suspect. The viscous stress tensor is assumed to be symmetric, and all 
variables are non-dimensional, see section 3.1 (24). 


For a large range of flight conditions, air may be assumed to be a perfect gas. This model 
condition of rigid rotating diatomic molecules may be represented by fixing y = 1.4, estimating 
static pressure from the equation of state, 

P - (y ♦ 1)p(* - ♦ ».] (2 6 > 
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Sutherland’s law is used to estimate the molecular dynamic coefficient of viscosity. 
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(2.8) 


Here T_ is the dimensional free-stream temperature, T « P/p, and thermal conductivity is estimated 
as k * T 0 71 . 


Equilibrium real gas calculations, which become important for the determination of body surface 
temperature and density at high mach numbers are an option available to the user of TEAM (17, 
25). These calculations are performed using curve fits developed by Srinivasan et al. [26, 27). The 
user must specify free-stream static pressure and density, which are used to estimate free-stream 
temperature, and free-stream specific enthalpy and specific energy, the ratio of which define y- 
Values for local static pressure, temperature, speed of sound, viscosity coefficient and thermal 
conductivity are then determined from the aforementioned curve fits using the estimated values of 
local specific energy and density. 

For viscous computations, the turbulence model used in TEAM is the well known Baldwin-Lomax 
Turbulence Model (BLTM) (28), or one can use laminar flow, lnviscid computations are carried 
out by solving the Euler equations. These will not be written explicitly here as they are derived 
by simply setting the right hand side of equation (2.1 ) equal to zero. The Euler equations may be 
solved with either perfect or real gas computation models as previously described. 

Kumerical dissipation needs to be included in the TEAM code primarily for two reasons, 1) 
stability of the solution process and 2) shock capturing [24). Stability must be numerically 
enhanced, as inviscid or high Reynolds number calculations have little or no physical dissipative 
phenomena. Without this artificial dissipation, the solution may become "saw-toothed" with 
alternating signs at neighboring cells. Often solutions of the RANS equations also require numerical 
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dissipation when the physical (viscous) dissipation is not adequate. 

fh?n^ a r P ha«H 8 J S performed automatically when solving the RANS equations. Euler solutions, on 

hw nf h ?L h rmi d C0 " tain the means 10 enforce an entropy condition as required bv the second 
law o f thermodynamics, hence providing a solution which is not phvsicallv realistic. The addition 
of dissipative terms which imitate the physics inside a shock wave circumvents this error TEAM 
«uriv de n S K a Ch0,Ce ° f ^ adapt,ve and characteristic based dissipation schemes. As thists a 
acetate Si 0 !'”’ 3 Chara . C,er,SllC ' based schcme is em P‘oyed. An upwind second order (USO) 
strong , °r n a ! Symme,r,C TVD Emulation (29. 30] allows for the capturing of 

;‘™ g ^ocks encountered at hypersonic speeds. The price of this improved shock capturing is an 

wssibk :^rr 0f H ar,thm r iC ° pera,ions ’ 35 spared to adaptive schemes. It is also not 
hvreclhzZ ,af v ‘ h p e condition of constant total enthalpy for a steady state Euler solution, which 

* P ST d b> h ,ncons,stency of ,he steady-state mass and energy conservation equations each 
of whose dissipative terms are constructed independently. q eacn 

TEAM requires the user to specify a boundary condition at all grid edges and interfaces Thic 

vatue^oTthe^eDe 1 !! dant^v L° * %h ° St Ce " S " bey0nd ,he arid boundaries and assjn a 

e of the dependant variable in the image cell which, when the fluxes of the cells (boundarv and 

aHL C Ulf fa V ce ra8ed ’ *"** Pr0Per boundary condition dependant evaluation of the flux vector 

Hl r n/r!nl^ Undary conditions are s P ecified at boundaries where the flow is incoming or outgoing 

diCt3,eS tha ‘, a " ° f thC f,0W quantit * es «»* image dll at St 
ooundanes be set to their free-stream values. At the outflow boundary all image cell Quantity 

“ ” W , the ' r T n<lary “ !l valu “- Th «' ■>« <te«rmined b he di “c oTof 

characteristics at the corresponding boundaries. uirection ot the 

mmeNav? ^"r dary co " dl " on5 are Prescribed differently for the Euler equations as compared 

l°soUd ^Tnd«.f VTam 0 " 5 L" !dd n ° W r ' 0JirB lhe "P-"“rmal-riow condition be satisfied 
at solid boundaries TEAM provides a choice of three methods to satisfy this condition That used 

for this research ,s the simplest and most robust of the three Surface pressureis set eoual o me 

,h ' "!! fiC ' " ,h ' ° n ' y Va " able “ comribu,e •» ">e momentum 

owcondmon A surface Z a ' *!■ ” “ r0 *' <»' C ' M face to preserve the no-normal- 

nnHYnn TI surface boundary condition for solution of the RANS equations i< the no-sliD 

ormomentum'I C h ,S r P ° Sed by SeUing lhe image - ce " of the cartesian'components 

l“Xe Thes^ m ,,e8a, :h e a 0 r ,he b0l,nda ' y -«" *»"-«• insuring ,he momentum be zero a 
f0 , “ , 1, c™, face m estimating surface pressure as in .he inviscid cases was employed 

he caMof.hela,,. ! """’"T" may b ' prescribed, or an adiabalic condition imposed. In 

va*u^^°he^magr- r ceM a mnt^ani a re temPtra,,Jre erad '"’' “ '"’' >0!ed » lk ' » «■•»» -he 

tolhtimaae “~n m “ S ' b ' !p *5' fied as ■ fluid ‘ conditions in the boundary condiiion dalasei. Values 

bicn‘ cui ° f iha b ^ da ^' iia — - 

Plane 11 ',? across' 'an V°? by m 'r'0'"'8 the flow field across she 

P , " ’ eg acr ®f an x_z P |ane of symmetry, the Y-component cf momentum changes sign while 

r - h ' if f bb “ ada ' y -«" ^">erppr.s. F„| me 2°. rn rue'll s m ri a 1 ca se s 

ed in this research effort, two planes of symmetry were specified with one cell between them. 
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(2.9) 


The semi-discrete approximation which is to be integrated in time is as follows. 

- <? v -£>*0 
at 

Here Q is the convective flux, Q v is the viscous flux, D is the dissipation and 0 represents the 
volume. Since the volume n is independent of time, equation (2.9) may be rewritten as, 

~~ * R(w) ' 0 (2.10) 

at 


where R is the residual defined as, 

- ^(Q c ( h>) - Q v (w) - D(w)) (2.11) 

Thus is defined a system of ordinary differential equations which may be solved by a variety of 
time marching schemes. Time accuracy is not important here as a computationally efficient steady 
state solution is the goal. TEAM uses an explicit multistage time-stepping scheme. This scheme 
allows relatively large time-steps and is easily vectorizable to exploit the capabilities of modern 
supercomputers. This m-stage hybrid scheme can be represented as follows. 

* M '* 

*,(1) m H .(0) _ 

*(2) = H ,(o> _ 


* w><°> - o B .jAr*R ( "‘ l) 

= h,<0) . 

M ,*«l _ H ,(») 

Where At* = CFL &t, the Courant number, CFL is a user specified parameter which scales the time 
step At. A pseudo time stepping, or spatially varying time step substantiallv reduces the number 
of time steps to convergence. This involves using a local time step for each cell, rather than a 
globally minimum time step. One consequence of pseudo time stepping is that the solution is no 
longer time accurate, ie. the solution is meaningless until convergence is reached. Viscous 
computations often require a much smaller time step than inviscid calculations because they require 
finer, highly clustered grids. TEAM allows the user to choose between three options for selecting 
the time step. 1) An inviscid time step in conjunction with two evaluations of numerical and 
viscous dissipation, 2) application of the modified Crocco’s scheme to scale the inviscid time step 
to satisfy a viscous stability limit, 3) use of a formulation proposed by Tannehill et al. [31], which 
estimates the time step to automatically satisfy the viscous stability criteria. In this research 
project, method 1) for inviscid calculations, and method 3) for viscous calculations has been used. 

Aerodynamic forces and moments on a body are determined by integrating the normal and 
tangential stresses. Shear stresses are, of course, absent for an inviscid computation. Denoting 
normal forces by superscript N, and shear forces by superscript S, force vectors and force 
coefficients in the body-fixed coordinate system are as follows. 
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Here p is static pressure, and 7 is the dot product of the stress tensor and the unit normal vector 

ft. 

The moments and moment coefficients about a point with position vector 7 are given by: 

J?>J l?* + A? 5 and 
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(32) 


Here A R is the reference area, and c R denotes the reference chord used in defining the force and 
moment coefficients. Coefficients of lift, drag and side-force, and the coefficients of the pitching, 
rolling and yawing moments are obtained in the wind-axis frame through the use of a 
transformation matrix. 



-S/Ate 

0 

COSa 

c n 

x 

COSpCOSo 

s/np 

COSpSWa 

c rr 


-SlNfiCOSa 

COSp 

- SIN $ SIN a 

l C «l 


CJ 


-sinP cosa 

cosP -sinj) sina 




r 

cocP coca 

sinp cocpsina 


C*r 

C ". 


-stna 

0 cosa 


C*i\ 


The expressions for aerodynamic parameters are for steady flow only. 
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